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We demonstrate the non-universal behavior of finite size scaling in (1+1) dimension of a nonlinear 
discrete growth model involving extended particles in generalized point of view. In particular, 
we show the violation of the universal nature of the scaling function corresponding to the height 
fluctuation in (1+1) dimension. The 2nd order moment of the height fluctuation shows three 
distinct crossover regions separated by two crossover time scales namely, txi and tx2. Each regime 
has different scaling property. The overall scaling behavior is postulated with a new scaling relation 
represented as the linear sum of two scaling functions valid for each scaling regime. Besides, we 
notice the dependence of the roughness exponents on the finite size of the system. The roughness 
exponents corresponding to the rough surface is compared with the growth rate or the velocity of 
the surface. 

PACS numbers; 68.35.Ct, 64.60.al, 61.43.Hv, 05.40.-a 



I. INTRODUCTION 

Theoretical and experimental study of morphology of the growing surfaces and interfaces, which are generated from 
various different growth processes, is a major challenge in recent years P, 01 @1- Several discrete models have been 
proposed to study the morphology of the rough surfaces formed in different growth processes such as imbibition in 
porous media 0, H, 0], thin film growth in molecular beam epitaxy (MBE) and growth of bacterial colonies 
These discrete models are based on simple stochastic growth rules, such as aggregation and diffusion. 

The rough surface evolved from a non-equilibrium process can be characterized from the study of the moments of 
the average of height fluctuations, called surface width, defined as 



W^iL,t) 



L 

i=l 



1/7 



(1) 



Where L is the system size, {h{t)) ^ X^i^i h{i,t) is the average height over different sites at time t and 7 is the 
order of the moment of the height fluctuation. 

The dynamic formation of the rough surfaces can be realized with the numerical simulation of discrete models. The 
dynamic and saturated behaviors of the rough surface are characterized with the power law nature of the moments 
of the height fluctuation. The short time height fluctuation is characterized by the exponent /3 known as dynamic 
one. The saturated behavior of the surface width determines the fractal dimension (represented in terms of roughness 
exponent a) of the rough surface. The overall scaling behavior of the rough surface is studied with the Family- Vicsek 
scaling ansatz (T^]. W2{L,t) being the 2nd order moment of the height fluctuation, one can write the Family- Vicsek 
phenomenological scaling ansatz as 

W2{L,t) = LV (^-1,^ (2) 

Where z = a/f3 and f{u) ^ [u <^1], f{u) ~ constant [w ;» 1]. This scaling relation is the central quantity of 
interest to study the morphology of the surface for any discrete growth process. The values of the exponents a and z 
uniquely determine the universality classes of the kinetic roughening process. Linear discrete model such as random 
deposition with surface relaxation (RDSR) [ll[ and nonlinear discrete models like ballistic deposition (BD) model [T^ . 
Eden growth (ED) model [l^, restricted solid-on-solid (RSOS) model and body-centered solid-on-solid (BCSOS) 
model [1^ have been proposed to study the kinetic roughening of the growth processes which belong to different 
universality classes. A competitive growth model consisting of two kinds of particles viz. BD (with probability 1—p) 
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and RDSR (with probability p) was reported . Two system dependent crossover time scales tc and r were found 
to characterize the rough surface. In this model, scaling properties follow RDSR model for < <C while for t 3> tc BD 
scaling properties are dominant over the concerned regime. The nonlinear coupling (A) is scaled with the abundance 
of the particles a,s X ^ p"^. The overall scaling was found to be 

W^2(L,t) - t'^«"^« i<tc 

W2{L, t) ~ (Ci + C2p"^W t » r (3) 

The suffixes correspond the respective models. Ci and C2 are constants, a is same for both the models. 

Several continuum models have also been prescribed to characterize the rough surfaces formed due to many natural 
growth processes ■ Continuum linear growth model described by Edward- Wilkinson (EW) equation [l^ belongs 
to the same universality class that of the discrete RDSR model. The values of a and /? in (1+1) dimension for this 
universality class are given as a ~ 1/2 and (3 = 1/4. The nonlinear Kardar-Parisi-Zhang (KPZ) equation [l^ defines 
the universality class which includes the discrete growth models, such as BD model, ED model, RSOS model and 
BCSOS model. In (1 + 1) dimension the scaling exponents for this class are a = 1/2 and (3 — 1/3 with the scaling 
identity 



a + a/f3 = 2 (4) 

This scaling identity follows due to Galilean invariance of the interface 12011 . However, nonlocality in the form of 
long range interaction in KPZ equation can modify the above identity [2l]. |22|. [23t . 

Instead of strong agreement between the theoretical and numerical predictions of a and j3 in the same universality 
class, the results do not match with several experimental findings. In (1+1) dimension, experiments on immiscible 
fluid displacement show the rou ghn ess exponent a lying between 0.73 and 0.89 ^24,, [2^, [S^] . Besides, experiment 
on growth of bacterial colonies 12811 yields a ~ 0.81. A possible explanation of the above mentioned experimental 
results was proposed by Zhang [29[ . He suggested a model with the noise amplitude having power law distribution as 
P{ri{'r,t)) ~ r]~^^^'^\ where r]{r,t) is delta-correlated noise. 

During recent decades, the growth of organic crystals on inorganic substrates draws great attention towards the 
application in electronic devices [s^l • Thin film growth of several organic crystals on different substrates (3l] . [s^ . 
have been studied experimentally. One diffusion-limited-aggregation (DLA) type of discrete growth model [3J] has 
been proposed to study the kinetic roughening of the growth of pentacene film on Si02 substrate. The growth of 
organic crystals involves particles of large sizes, viz, the molecular formula of pentacene, copper phthalocyanin and 
3,4,9,10-perylenetetracarboxylic dianhydride are C22H14, C32HiqNsCu and C24HsOq respectively. 

Motivated by the above various growth mechanisms which involve small and extended particles with different types 
of relaxation rules, we have proposed a discrete growth model in (1+1) dimension. This model may involve particles 
of different sizes in a single growth process. In this model, the particle size plays a major role in the kinetic roughening 
of the surface. This model shows a morphological transition from multifractal to unifractal regime beyond a system 
dependent characteristic length scale [35 1. The morphology of the surface depends on the dominating particles in the 
growth process. The finite size scaling relation is not of the universal type defined in equation not even satisfying 
the equation A new scaling relation is proposed to characterize the rough surface. The values of the scaling 
exponents a and (3 are not universal for this model. It has a finite size dependence with a particular scaling form. 

The paper is organized as follows. In Sec II we describe the model with the simulation parameters. In Sec HI, we 
propose the scaling law characterizing the rough surface. The effect of finite size on the scaling exponents is discussed 
in section IV. Porosity, a parameter that represents the bulk nature of the system, is introduced in a new scaling form 
in section V. Finally, in Sec VI the conclusions are drawn from the numerical results described in previous sections. 



II. THE MODEL AND ITS SIMULATION 



The model introduced here contains essentially the aggregation and diffusion mechanisms with a new idea of 
extended particles from a generalized point of view. The concept of these participating particles is introduced in such 
a way so that the particles can be of various sizes for a single growth process. The diffusion mechanism considered 
here is same for all types of particles. In spite of the same diffusion mechanism for all particles, however the final 
morphology of the surface is not the same when particles of different sizes are taking part in the growth process 
separately. A detailed description of the model is given below. 

The substrate lattice is taken of size L with periodic boundary condition. The participating particles in this model 
are considered as different multiple matrix sequence in terms of the smallest unit of the substrate lattice. In this way, 
if the substrate lattice is considered as a 1 x L matrix such that the smallest unit of this lattice is a 1 x 1 matrix , 
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then the involving particles may be of sizes as 1 x 

1 (□), 1 X 2 (□ ^ HO), 2 X 1 (U = □ )and so on. From 
this generalized point of view, we call the particles as 'extended particle'. The units of 1 x 1 particles at the ends 
of each extended particles which are facing towards the substrate lattice will be called 'extreme cells'. In this sense, 

the particles 1 x 1, 2 x 1, 3 x 1, have only one extreme cell while the particles 1 x 2, 2 x 2, 1 x 3, have two 

extreme cells (see Figll] extreme cells are shown as shaded boxes). 
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FIG. 1: Different types of particles participating the growth process. (1) is the 1x1 particle, (2) is the 2x1 particle, (3) is 
the 1x2 particle and so on. The shaded boxes of each particle are the extreme cells as defined in the text. 

When these above mentioned particles form the rough surface, a 'stable position' should be maintained throughout 
the surface for each of the involving particles. The stable position refers to a condition when at least one point from 
each of the extreme cells for each extended particles get contact with the columns from the substrate lattice separately. 
In other words, the extended particles will attach with the substrate lattice if at least one corner from each extreme 
cell is shared with the corners of the columns from the substrate lattice. In this circumstance, the selection of the cell 
from the two extreme cells by the substrate lattice site is random. The overall aggregation and diffusion mechanism 
of the different types of particles is described as follows: one of the extreme cells of each of the extended particles is 
chosen randomly by a substrate site, then it (the extended particle) slide according to the lower height profile of the 
nearest-neighbor columns with a stable position. The process continues till the stable position is reached. Naturally, 
when the the number of extreme cells is one, then the selection of this extreme cell by the substrate lattice site is 
completely deterministic. 

It is obvious that the aggregation and diffusion mechanism of the 1x1 particles follow the RDSR mechanism [T]| . 
The 2x1 particles also follows the RDSR mechanism except the height increment of the column at which it finally 
sticks is of two units rather one as that of 1 x 1 particles. Similar rules are also valid for the relaxation mechanism of the 

particles of sizes 3x1,4x1, But the case is quite different for the 1x2,2x2, 1x3, particles. The mechanism 

of the diffusion of these particles is not that of RDSR type. The growth rules are straightforward generalization of the 
RDSR model for extended particles. It seems that the diffusion of each of such type of particles will cause different 
type of surface morphology. For clarification, we show in Figl2]and Fig[3] schematically the aggregation and diffusion 
mechanism for the 1x2 and 1x3 particles respectively. In both of these two figures, the one of the extreme cells (I 
and II) corresponding to the extended particles is chosen by a site (A) of the substrate lattice randomly. 

With these rules, we have developed a model in (1 + 1) dimension with the participating particles having different 
sizes. The morphology of the surface should be determined by the particles dominating the growth process. During 
the growth process bulk defects are allowed to form. But the growth rules have been set in such a way that the 
voids formed in the system are closed, so that possible overhangs should be avoided. The lateral growth property 
breaks the up-down symmetry, resulting the presence of non-linearity due to the local slope (Vft-(r, t)) fluctuation in 
the continuum description of the model. It was previously reported [1^ that the participation of the 2x1 particles 
in the growth process does not significantly change the morphology of the surface generated from RDSR mechanism. 
Moreover, the 1x2 particles create voids in the system during surface growth, so KPZ type of growth characteristics 
may occur. Also the up-down symmetry is broken due to the involvement of the 1x2 particles in the growth process. 
So, the scaling properties of the kinetic roughening of the surface are expected to be modified by the non-linearity 
introduced by 1 x 2 particles. Due to the above mentioned reasons we have simulated the surface and analyzed its 
kinetic roughening for the surface formed due to the 1x2 particles only. The morphology of the surface formed due 
to the deposition of 1 x 2 particles for the system size L = 200 at time 2 x lO'^ is shown in Fig[4]with four equal time 
zones. The increment of roughness of the surface with time is evident from the figure. 

The model has been simulated in the lattices of lengths L = 50 x 2" (n = to n = 6). Deposition time is taken 
as 10^ < t < 10'' depending upon the lattice size. The internal structure of the surface is characterized by 'intrinsic 
width' [3^. The probable origin of intrinsic width is the voids, overhangs and large local slopes. To minimize the 
effects from intrinsic width, we have incorporated the noise reduction technique [s^. The noise reduction parameter is 
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FIG. 2: The aggregation and diffusion mechanism for 1x2 particles. Different possible cases are shown. Here one of the 
extreme cell (I) is chosen by a site (A) of the substrate lattice randomly. The case will be equivalent when the other extreme 
cell (II) will be chosen by the site (A) of the substrate lattice. 




FIG. 3: The deposition and relaxation mechanism for the 1x3 particles. Different possible cases are shown. A random lattice 
site (A) chooses one of the extreme cells (I) of the extended 1x3 particle randomly. The choice of another extreme cell (II) 
give a statistically same type aggregation and diffusion mechanism. 



the number of attempts per site for the actual aggregation process. In this model, the noise reduction parameter (to) 
was set fixed as to = 10. With this value of m the surface morphology shows a stable scaling behavior with repeated 
independent simulations. Uniformly distributed uncorrelated noise has been taken. Depending upon the system size, 
the results were averaged over 100 to 10 independent runs. Simulations were done on an IBM Server PC with two 
64-bit quad-core P0WER5+ processors. 
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FIG. 4: The rough surface formed due to the deposition of 1 x 2 particles. The total number of particles involved here, have 
been isolated into four set with different color shading showing the change in roughness. 



III. NON-UNIVERSAL SCALING OF KINETIC ROUGHENING OF THE SURFACES 

The kinetic roughening is characterized by the scahng exponents corresponding to the moments of height fluctuation 
defined in equation ([T]) . The log-log plot of the 2nd order moment of the height fluctuation for the values of lattice 
size 50 < L < 3200 is shown in FigO It is clear from the figure that three distinct scaling regimes exist, which are 
separated by two time scales namely t^i and tx2- The situation occurred around the crossover time scale t^i should 
not be confused with the situation occurred when the finite size scaling is affected by the intrinsic width of the system. 
Because intrinsic width is a system size independent effect (STj . The crossover between the intrinsic width affected 
regime with the unaffected regime is independent of the system size. But from Figl5]it is clear that the crossover txi 
is system size dependent. The scaling form of the time scale txi will be shown latter. 
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FIG. 5: Log-log plot of W2{L,t), the 2nd order moment of the height fluctuation as defined from equation ([T]), versus t. The 
dotted lines are the power law fitted according to the equation ([5]) in three distinct regimes. 

To demonstrate the sealing behavior of the rough surface evolved from the present discrete model, the time evolution 
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FIG. 6: Representation of different scaling regimes of the 2nd order height fluctuation for a particular system size L=400 with 
the respective power law fitting. This fit shows a = 0.698 ± 0.002, /3i = 0.22 ± 0.02 and /32 = 1.21 ± 0.04. 



of the height fluctuation corresponding to a particular system size {L = 400) has been shown in FiglH The nature 
of the plot in FiglBlis quite similar to that of the plot for the competitive growth model consisting of RDSR and BD 
[l6j . But the situation is different here. Only one kind of particles arc taking part in the growth process. So. no 
other time scale can appear depending upon the abundances of the different kinds of particles participating in the 
growth process as in the discrete model described in ref [iGj . Thus, the scaling behavior for the present model can 
not be defined by the equation ([3]) . Guided by the nature of the scaling of a growth process [ij , three power laws with 
different exponents have been fitted in three different scaling regimes as follows 
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Also, in FigOthe power law fitting for different system sizes has been shown by dotted lines. The scaling behavior 
shown in equation ([5]) indicates that, though the surface evolved with a definite mechanism, two different growth 
processes belonging to different universality classes are governing the overall kinetic roughening process. Such a scaling 
behavior of a growth model has not been noticed previously in literature. The system contains three characteristic 
time scales <xi, ^ and tx2 (see FiglS]). The time scale r denotes the saturation time for a growth process belonging 
to a certain universality class with the growth exponents a and (3i. Another growth process belonging to different 
universality class seems to occur beyond t ^ t^i with the growth exponents a and ^2- From FiglHl it is seen that 
the two independent growth processes having saturation time scales r and ix2 actually determine the crossover time 
scale txi- So, the crossover time scale t^i should depend on the time scales r and ix2- The crucial time scale t^i is 
the crossover between these two growth processes. The saturated rough surface is characterized by the latter growth 
process dominating beyond txi- The time scales corresponding to the two distinguished growth processes should 
behave as 

ix2 ~ L^' (6) 

Also from the nature of the height fluctuation (see FiglS]) it can be argued that the crossover time scale txi can be 
scaled with the system size L as 

txi-L^ (7) 

Where zi = a/Pi and Z2 ^ a/ (32 are the two dynamic exponents corresponding to the two growth processes occurring 
in two different time regimes, x is a different scaling exponent. Since the crossover time scale txi is dependent on the 
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other two independent time scales t and fxi so the scahng exponents x depends naturally on the two independent 
exponents zi and Z2- 

From the scaling behavior, we argue the simplest possibility that the kinetic roughening of the surface is occurring 
by the two growth processes in two different time regimes (viz t ^ txi and t ^ txi) independently. With this 
consideration, we propose that the overall morphology of the surface is governed by the scaling relation which can 
be represented as the linear sum of two scaling functions corresponding to each of the growth process dominating in 
different regimes, leading to the saturated rough surface having unique roughness exponent. So, mathematically the 
proposed scaling relation can be represented as 



where the various scaling functions are defined by 



(8) 



fi (u^ ) ^ constant u-i 3> 

f2{u2) ^ constant U2 <S — _ 

X 

f2{u2) ~ constant U2 ^ 1 (9) 

To observe the appropriate scaling and the crossover, we proceed as follows. According to the scaling relation ([5]), 
in the time regime t <S txi i-e, for t/L^^ <C 1/L~^~^ the scaling functions /i and /2 will be 



/: 
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Thus from equation ([5]) the scaling relation becomes 



/2 ( — — I ^ constant 



t \'^^ t 1 
W^2(i,t)/i" ~ ( ^ when — <- (10) 



For the time regime <xi -C t -C tx2 i-e, for ^ <^ t/L^'^ <C 1 the scaling functions /i and /2 will looks like 

fi I I ^ constant 

t\ f t 



So according to the equation ([8]) the scaling relation reduces to 



t \^-' . 1 t 



W2{L,t)/L'^^^—j when — (11) 

The time regime t ^ tx2 i-e, when t/L^'^ ^ 1 the scaling fimctions /i and /2 behave like 

/i (^^T'^ ^ constant 
t 



f2 ( -r^ ] ^ constant 
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In this time regime the scahng relation turns out as 



W2{L,t)/L" constant when 7^ > 1 (12) 

J_/ 2 



The above scaling relations in equations (jlOp . ([TT|) and (|12p satisfy the observation in equation ([5]). The complete 
data collapse in the log-log plot of W2(i,t)/i" in different time regimes as shown the Figl7]also confirm the above 
scaling relation. 

Thus, for this present discrete model, the kinetic roughening of the surface can not be scaled with a unique scaling 
function having two independent scaling exponents. It thus deviates from the universal nature of the scaling of the 
kinetic roughening of the surfaces. 

The visualization of the scaling relation defined in equation ([8]) with the scaling functions given in equation ([9]), is 
shown in Figl7]with the data collapse in two different time regimes t ^txi and t 3> ixi- The scaling functions have 
been plotted for different set of L, a, /3i and P2 values, because of the strong finite size effect on a and weak finite 
size effect on (3i and (32- Though the x-coordinatc for each plot arc different (one is t/L^^ and the other t/L^^), we 
plot in the same graph to show the existence of two distinguished scaling regimes t ^ txi and t ^ txi- The large 
gap between the two scaling regimes, t txi and t ^ txi, is due to the large difference of two dynamic exponents 
zi and 22- Also, long crossover region around t ^ txi is one of the reason for such a large gap. The data collapse for 
scaling functions defined in equation ([9]) with the exponents defined in equation ([5]) for the respective time regimes 
defined in equations ^ and (O, confirms the scaling relation defined by the equation ([5]). 




FIG. 7: The scaling plot in log-log scale shows that W-ziL, t)/L°' is plotted against t/L^^ for t < L'^^ and WziL, t) /L" is plotted 
against t/L'^ for t ^ . As predicted from the scaling relation given by the equation ((S)), two distinct scaling functions given 
by equation Q are seen to exist prominently. 

The scaling behavior shown in FigO can be explained physically in the following way. In our model, the particles 
diffuse along the surface in search of a stable position. During the growth process, there is a interplay of diffusion 
with the finding the stable position for each particle. Initially, since the total number of involved particles in the 
growth process is small, each particle is allowed to diffuse throughout the surface. So, within this time, diffusion of 
particles is dominated over the process of finding the stable position. After the relaxation, when a 1 x 2 particle 
sticks to the substrate lattice, it implies that the constituent two 1x1 particles stick with the same height. Thus, 
from the 1x1 particle point of view the height-height fluctuation also decreases. Due to the above two reasons the 
overall height-height fluctuation of the surface during this time period will be very small, resulting a small value of 
the dynamic exponent However with the increase of time, more and more particles do take part in the growth 
process. Due to the relaxation rules there is a fair probability of getting a stable position for each particle within 
the short-range sites. Effectively, the particles are restricted to the connected sites of the lattice site on which they 
were deposited. In other words, the particles are now localized almost stopping the long-range relaxation. Thus, an 
instability occurs due to the piling of particles on the upper terraces and restriction of the relaxation on the lower 
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terraces. This increases the height- height fluctuation, resulting a large vahie of the dynamic exponent (/3). This effect 
is quite similar to the Ehrlich-Schwoebel (ES) effect [3^, HO, [4l|, [13, IS^ which arises due the presence of a ES barrier. 
This effect induces an instability by hindering step-edge atoms on upper terraces from going down to lower terraces 
in the MBE type of growth processes. A discrete solid-on-solid model [44| of epitaxial growth without bulk defect was 
proposed, which takes into account the ES effect with the introduction of a parameter (Rmr) called the incorporation 
radius as a ES barrier. Another reasonably realistic bulk defect induced discrete model [4^ of epitaxial growth in 
(1 -|- 1) dimension was presented in such a way so that the kinetic roughening is controlled by the interplay of the 
mound instability with the KPZ roughening. In this model, the diffusion of the particles were partially controlled by 
the parameter which actually selects the direction of diffusion with a probability exp(— i?). In both of the above 
models the instability occurs due to the presence of a step edge barrier, in former case it is infinite while in the latter 
case it is finite. However, in our model, we do not put any step edge barrier explicitly, the instability occurred here 
is completely self organized. 

The complicated scaling behavior of the rough surface evolved from the present discrete model is represented with 
a linear sum of two independent scaling functions corresponding to different growth processes. Now we point out 
the unusual behavior of the scaling function around t <xi- The crossover around txi leads to a morphological 
'phase transition' from one universality class to another universality class. Such a morphological phase transition 
with a unique growth mechanism has not been observed previously for any discrete growth model. In this context, 
we may also mention that another morphological linear-nonlinear 'phase transition' is seen to occur around a critical 
probability of deposition of the 1x2 particles, beyond a characteristic length scale, for a competitive growth model 
involving particles of sizes 1x1,2x1 and 1x2 [sl] . Below we would like to study the variation of the roughness 
exponents with the system size L. 



IV. FINITE SIZE EFFECT ON ROUGHNESS EXPONENT 

The earlier work by Krug and Meakin [s^l had shown that for a nonlinear KPZ growth model the roughness 
exponents were affected by the finite size dependence of the steady state growth rate of the system. The demonstration 
of the finite size effect on the growth rate of the evolution of the surface is shown first. The growth rate of the surface 
is defined as V{L,t) = According to the suggestion by Krug and Meakin [s^, the steady state growth rate is 

scaled with the system size L as 

T/(i,t ^ oo) = ^ oo) - Ai-'" (13) 
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FIG. 8: Fitting of steady state velocity in the asymptotic limit versus L " with v = 0.61 ± 0.05. 

The scaling behavior of the velocity V in equation (fT3|) is shown in FigEl The fitting shows v — 0.61 ± 0.05, with 
asymptotic limit of the steady state velocity is V{L oo) = 1.083 ± 0.02. It was predicted [111 that the roughness 
exponent for large enough system for growth models belongs to the KPZ universality class, would be related with v 
as 

a = \-vl2 (14) 
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Strong finite size effect on a is also observed in the present model, 
dependence on a has been suggested [i^ earlier. 



The following scaling relation for the finite size 



a{L) = a{L ^ oo) + TL 



-5 



(15) 



In Figini the a values are fitted according to the above scaling relation given by equation (fT5)l with 6 = 0.57 ±0.03. 
For asymptotically infinite system the roughness exponent a{L — > oo) = 0.794 ± 0.005 which is comparable with the 
prediction in equation (fT4|) . The relation between the exponents shown in equation (fT4|) is based on the realization 
that the Family- Vicsck scaling relation (equation ^) is satisfied. In the hydrodynamic limit, the present model 
satisfies the scaling relation given in equation ^ rather than in equation ([2]). So, it can be argued that the correction 
due to the finite size dependence of the roughness exponent will be not like that of the model which follows KPZ type 
of growth. We have compared these two models because, both of these models have inclination dependent current 
due to the lateral growth property, which is the source of KPZ like behavior in the grow th process. Growth models 
with power law distributed noise events P{ri{r,t)) ~ show such values of a j47l. [ist with different values of 

fj,. However, in our present model, the noise distribution is (5-correlated with uniform amplitude. Rare events are not 
occurring here. With the time evolution of the system, a multifractal behavior is also seen to occur in the present 
model [3a | , similar in nature as that of the rare event dominated growth model (49j . 




FIG. 9: Finite size dependence of the roughness exponent a. The fitting of a with L shows that S = 0.57 ± 0.03. 

The values of the two dynamic exponents /3i and /?2 remains steady for large system size. The values of this two 
dynamic exponents are found to be /?i = 0.22 ± 0.03 and /?2 = 1.22 ± 0.01 for large enough systems. Such a large 
value of the dynamic exponent /J2 still has not been observed in literature. Due to strong finite size dependence of a 
and weak finite size effect on the values of (3i and P2, the crossover time scale txi dependent on system size L with a 
exponent x given in equation ([7]). Moreover, a + zi 7^ 2 and a + Z2 7^ 2; this immediately points out the breakdown 
of the Galilean invariance of the system. It is to be noted that the value of dynamic exponent /3i is very close to that 
of the EW model with /3i = i. Above the system dependent time scale < ~ txi the dynamic exponent /32 has a very 
high value. The transition of the dynamic exponent from low value to a higher value can be visualized by looking 
critically to the individual configuration of the surface and the bulk. As per the aggregation and diffusion rules, it 
appears that initially the surface moves compactly without having voids and that is why the dynamic exponent 
in that region is close to that of the RDSR model. In latter time, voids are incorporated into the system and an 
ES like instability is found to occur. This self organized ES effect triggers the rapid roughening of the surface. The 
dynamic exponent {(32) becomes high in that region. The experimental observation of rapid roughening involving 
extended particles in the growth of organic thin film (Diindenoperylene) was reported previously [50| ■ 



V. 



BULK PROPERTIES AND ITS SCALING 



To have a deep insight into the internal structure of the interface, the bulk properties of the system are of great 
interest. The diffusion mechanism for the present growth model is unique by its definition. Closed voids created due 
to such diffusion mechanism made the system porous in its own way. So, the bulk property will be different from 
the other porous systems, created from different kind of aggregation and diffusion mechanism. The bulk properties 
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of a system can be quantified with the definition of porosity. Since the number of closed voids can be determined 
accurately, we define, porosity P for this particular system in a similar way as defined in Ref [5]| . 



P = NjNt (16) 

Where — Number of voids and Nt — Number of voids + Number of particles deposited. A deposition process 
where the number of particles is conserved, is characterized by the particle flux J. As shown by Krug [S^], the deposit 
density p can be related with the growth velocity v as 



J/v 



(17) 



According to the equation pT]) the deposit density, which actually characterized by the quantity porosity, will be 
scaled with the system size as that of the velocity with the modification due to the particle flux dependence on the 
system size in the asymptotic limit. To see the asymptotic limit of the porous structure of the system at saturation 
with system size, we propose a scaling relation of the porosity as 



P{L, t ^ oo) = P{L ^ oo) + TL-'^ 
The above proposed scaling relation is new in the present literature. 



(18) 



0.12 



0.11 



I 0.10 
J 



J! * ■ 




, -0.72 



FIG. 10: Plot of porosity in the asymptotic limit with L ^ with the exponent 77 — 0.72 ± 0.02. 

From physical point view, porosity can be dcflned as the reverse of the deposit density. That is 

P^l-p (19) 

From equation (|17p for uniform flux J 

1-P = J/v 

P (20) 

Fig) 101 shows the behavior of the porosity at the asymptotic limit with different system sizes. It shows that asymptotic 
value of the porosity P{L 00) as 0.087 ± 0.005 and the coefficient F = 0.34 ± 0.013. The prediction from equation 
([20|) is well in agreement with the results shown in FigEJand FigfTOl 

VI. CONCLUSION 

The kinetic roughening of the surface created due to a nonlinear discrete growth model, is studied here. Several 
features, not previously observed corresponding to kinetic roughening, are observed in the present model in (1+1) 
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dimension. To summarize, we mention the following points systematically. The finite size scaling of the rough surface 
shows a different type scaling nature. Two distinguished time scales, corresponding to the height fluctuation, emerge 
in the system in (1+1) dimension. They separate three scaling regimes with different scaling exponents as well as 
scaling functions. To characterize this kinetic roughening, a new scaling relation is proposed, which is represented 
as the linear sum of two scaling functions valid for two distinguished scaling regimes. The existence of two scaling 
regimes with small and large values of dynamic exponent /3 is well explained with the occurrence of a self organized 
Ehrlich-Schwoebel (ES) like instability (caused due to the localization of the extended particles) which triggers the 
rapid roughening of the surface. Due to the finite size effect on the growth rate, the scaling exponents are also affected 
by the finite size of the system. The finite size effect on the roughness exponent is scaled with a scaling relation. The 
scaling exponent for this scaling relation is well compared with the prediction made by Krug and Meakin. The bulk 
nature of the system for different sizes is shown through a new scaling relation. 
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